Time-evolution of the Rule 150 cellular automaton activity from a Fibonacci iteration 
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The total activity of the single-seeded cellular rule 150 automaton does not follow a one-step 
iteration like other elementary cellular automata, but can be solved as a two-step vectorial, or 
string, iteration, which can be viewed as a generalization of Fibonacci iteration generating the time 
series from a sequence of vectors of increasing length. This allows to compute the total activity time 
series more efficiently than by simulating the whole spatio-temporal process, or even by using the 
closed expression. 

PACS numbers: 05.40.-a, 05.45.Df, 82.40.Np, 45.70.Qj 



I. INTRODUCTION 

Since the coining paper of Bak, Tang, and Wiesenfeld 
0, there has been considerable interest in the long-time 
behavior of cellular automata, especially for occurrence 
of long range correlations, and correspondingly for power 
spectra exhibiting a power law decay, which have become 
a paradigm for comi □lex dynamical systems in general 0. 
Cellular automata 0, are a widely studied class of dis- 
crete dynamical systems showing emergence of complex 
spatio-temporal patterns from a simple dynamical rule. 

ECA and sum signals — A cellular automaton con- 
sists of an infinite lattice of cells of e.g. two possible states 
(0,1), and a local deterministic update rule. At each dis- 



crete time step, a cell is updated x n 



,./ + ! 



according 



to the state within a local neighborhood. For Conway's 
Game of Life [f| the 3x3 Moore neighborhood on a two- 
dimensional lattice is used. A simpler, yet complex, class 
are the elementary cellular automata (EC A) 0, 3 , de- 
fined on a one-dimensional lattice, and the update rule 
depends on the next-neighbor sites and its own state one 
time step before: 
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where / (the rule) is determined by 8 bits being the out- 
put of the possible input bits 000, 001, Ill; this 8-bit 
number forms the "rule" number which enumerates the 
256 possible ECA rules. The power spectra of rule 90 
@ and some of the other ECA rules Q exhibit a 1/ f a 
decay. Rule 90 and rule 150 can be expressed also as 
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where r = defines rule 90, and r = 1 rule 150, respec- 
tively. In the context of catalytic processes, both can be 
interpreted as local self- limiting reaction processes 0, 0] • 
As models for chemical turbulence, similar dynamics with 
a continuous phase variable have been discussed in Q and 
0, including solitonic behavior, periodic, and turbulent 
states. 

Total activity. - In the chemical picture, the total 
reaction rate at a given time t corresponds to the total 
number of sites with x l m = 1, described by the sum signal 



(3) 



While for rule 150 there is no convenient solution of X(i) 
except a formal one |l5j | , and a fairly complicated expres- 
sion (see Sec.[nJ|, ^ is computationally quite costly to per- 
form the full spatiotemporal dynamical simulation, even 
if one is interested only in the time series. This paper 
gives an iterative solution of X(t) from a geometrical it- 
eration and investigates the relationship to the Fibonacci 
iteration. The block sums over < t < 2 N — 1 can even 
be expressed directly via Fibonacci numbers. 

Throughout this paper the pure pattern generated by 
a single 1 are considered on an infinite lattice. 

II. EXACT SOLUTION 

While in the Sierpinski (rule 90) case X(t) factorizes 
in a product of Xi(ai) for all "time spins" @, for rule 
150 it does not. But again, a "spin decomposition" of 
time t = X^jL^ 1 with aj G {0, 1} can be utilized as 
an efficient coordinate system for the time axis. 

Before turning to the geometric iteration, it should be 
mentioned that a closed expression in fact can be written 
down as follows. As pointed out by Wolfram Q, for 
rule 150 the "correlation" of the time spins comes into 
play, i.e. X(t) is exactly multiplicative for blocks of spins 
of value 1 which are separated by one or more zeroes. 
Then X(t) — Y[n=iX-( n )° n > where c„ is the multiplicity 
of blocks of length n. The series x( n ) should correctly 
read 0] to the iteration x( n ) = ^x( n — 1) — ( — 1)" f° r 
n < 1 and x(0) = 1. Obviously x(n) = X{2 N - 1) holds. 
The x{ n ) on t he other hand turn out as the most decaying 
frequencies in the spectrum of rule 90 (see Fig. 2 in Ref. 
0), and can be expressed in closed form by 



2' 
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where |_ J is the floor function. Defining <7_i := and 
<jn :— 0, one can formalize the spin-block counting as 

N—n n— 1 

N ^(l-^-OCl-^i+n) II ai + l 

X(t) = Hx(n)-° i= ° . (5) 

n=l 

With our expression this is a closed solution, and cor- 
responds to Xgo(t) = [j.Q 2 ffj in the Sierpinski case. 
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Due to the complicated time spin correlations, it how- 
ever looks quite unwieldy for analytical use, and even is 
numerically unfeasible []J| ■ 

III. ITERATIVE SOLUTION BY 
GENERALIZED HYPER-FIBONACCI SERIES 

In contrast to the rule 90 (Sierpinski) case, for rule 
150 the time evolution does not follow the same type 



J 



of initiator-generator mechanism as it is well known for 
fractal sets. However, it is possible to define a geometric 
or measure-theoretic iteration based on the last and 
the last but one iterate, see Fig.^ This corresponds to a 
difference equation with the r.h.s. depending on the last 
two time steps, and in fact, for the total activity within 
2 N time steps we will derive a difference equation later. 




FIG. 1: Left: Time-evolution of rule 150 for the first 64 time steps, started with a single seed. Right: Illustration of the 
replication rule. The whole system is symmetric with respect to the vertical axis. The whole triangle above is replicated once 
(left part of the triangle). The upper part is reproduced quadruplicate (right part of the triangle). 



According to the replication law (see Fig. QJ, the time series of the total activity X(t) follows the two-step iteration 



If we define 



the iteration reads 



X(l,l,cT n _ 2 ---CTo) = X(0,l,a n _ 2 ---o-o)+2X(0,0,CT Ti _ 2 ---<T ) 



X(1,0, (7 n _ 2 ■■■CTq) 



Z(l, er„-2 ' 
F(l,a„_ 2 . 



•ooj 



3X(0,0,a„_ 2 ---a ). 



X(.,0, 
X(m,l, 



Z(l, er„_2 • • • cr ) + 2Y(1, a n -2 ■■■o- ) 
+ 3F(l,cr„_2 • • - co), 

I 



(6) 
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or short 

K = 
K = 

and the concatenations 
Y n = 



3Y„_! 



(Y n _i, Z„_i) T 
Z« = (Z°,Z^) T 
complete the iteration 



and 



Y„J = 

Zn) = 



3 1 
21 1 



Y n _i 
Z n _i 

Y n _r 



(9) 

(10) 
(11) 

(12) 
(13) 



The initial vector is given by 



Y 



. Eqns. 



(I12I13|) can be collected together to the iteration 
Y„ 




(14) 



where the dimension of the vectors Y„, Z„ is 2 n , growing 
in the same way as for the Sierpinski iteration (see Fig- El) 
1 \ ( Y n _r 
2 • 1 ) ~ \ Z„_x 

1 \( Y„_! 

(-l)-l ) I Z„_ a 



eration ( Z„ I = 



or the Thue-Morse it- 
(both use O 



3 



and start with Zq = (1)). 




FIG. 2: Illustration of the replication rule for rule 90. 



or Lucas iteration 



Here 



These iterations look formally similar to the Fibonacci 
F n \_ f 1 1 \ f F n _i 
F„_! J \ 1 J \ F„_ 2 
(Fi,Fo) T = (1,0) T defines the Fibonacci series and 
(Fi,Fo) T = (1,2) T is the initial condition of the Lucas 
series. For the latter two, the length if the iterates is not 
growing. 

Eqns. |JBJ| or (|14f> . equivalently, together with the ini- 
tial condition (1,3) T generate X(t) iteratively for all t. 
Formally this iteration is analogeous to the Fibonacci it- 
eration, but acts on vectors of growing length within an 
infinite-dimensional vector space indexed by nonnegative 
integer values. 

This type of series should be distinguished from the 
(r-th) hyper-Fibonacci series H3, where ft = 2 1 - 1 for 
/ = 1, r + 1 and / ; = //_! + ••• + fi- r +i for I > r + 1. 
On the other hand, the terminus generalized Fibonacci 
series is widely used for the ordinary Fibonacci or Lucas 
iteration with two arbitrary start values fo and f\ , where 
fo = 0, fi = 1 defines the Fibonacci series and fo = 2, 
fi = 1 defines the Lucas series; and in fact both can 
be used as (nonorthogonal) basis vectors of the linear 
space of generalized Fibonacci series. The r = 1 hyper- 
Fibonacci series corresponds to a generalized Fibonacci 
series with /o = 1, /i = 2. Consequently, an iteration of 
the algebraic structure of Eq. (|14fl could be denoted as a 
generalized hyper-Fibonacci series. 

Another observation is the partial self-similarity relation 

X(o n ■ ■ ■ cr 3 , 0, cri, Co) = X(a-L,a ) ■ X(a n ■ ■ -a 3 ) (15) 

(leading zeroes omitted in notation), i.e. the sequence 
generated by every second four-block {a 2 = 0) factorizes 
into the first block (1, 3, 3, 5) and the whole sequence it- 
self. A closed expression for X(a n ■ ■ ■ 03, 1, ai, <jq) is how- 
ever not known yet. 

The first values of X(t) = X(t 1 +t 2 ) are listed in Tab.|l| 
(see Fig. 01. 



IV. ITERATIVE SOLUTION VIA A FORMAL 
LANGUAGE REPLICATION RULE 



One can proceed similar to the Sierpinski (rule 90) 
case, where X(t) can be generated from an initiator 
a = (1) by the iterative replication (a) — ► (a, 2 • a). 
For the Thue-Morse se quen ce, the iteration reads (a) — > 
(a, (—1) • a), see [l2llT^.lT^ . However, for rule 150 again 
the Fibonacci-like structure comes into play, i.e. each iter- 
ate depends on two preceders for which one has to require 
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TABLE I: The total activity X(t) for the first 256 time steps. 




FIG. 3: Rule 150: Plot of X(t) for the first 256 time steps. 

len(a) = len(fe). Then the iteration is defined by 

(a, b) -> (a, b, 3a, 2a + b). (16) 
This is equivalent to Eq. (|14f> . 



V. BLOCK-SUMS AND THE FIBONACCI 
SERIES 

Following the same geometrical argument as for the 

row sums, the sum S n = X)?=o 1 i s given by the 
iteration S n — S n -i — S n -i + 4S' n _2, or 

S n = 2S n - 1 +4S n -2, (17) 

and the first elements of the series are listed in Tab. [H] 
The matrix of the iteration in time-delayed coordinates 



Sn-i 




(18) 



has the eigenvalues — 1 ±v5, indicating that it dif- 
fers from the Fibonacci iteration matrix by an additional 
expansion factor of 2, i.e. by a suitable transformation 
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TABLE II: The block sums S n for the first 18 time steps. 



interval [0, 2^ — 1]. — The first values of this series are 
JV ,iVi,...= 1,3,4,7,11,18,.... 

An immediate side result is that S n defines the to- 
tal volume of sites in (n, t) space having the value 1, and 
this volume scales for n — > oo with the largest eigenvalue. 
Thus, if the time-doubling iteration is interpreted as gen- 
eration rule of the resulting self-similar fractal (rescaled 
to the unit interval), its Hausdorff-Besicovic dimension is 
given by (1 W5)/2. 



it relates to the Fibonacci numbers in usual convention, 



VI. CONCLUSIONS 



S n 



F 



n+2 



1 

V5 



1 + V5 N 



1 - V5 N 



(20) 
(21) 



the latter following from Binet's formula. A "blockwise 
normalization", or detrended signal, as used in |(|, can 
be acheived by subtracting 



N„ 



(S n — S n -i)/2 r 



2F n -2 — F n -2 



(22) 



fromX(f) for all 2"- 1 <t < 2 n -l. Then X(t) - N [log2 tJ 
defines a signal with a mean vanishing within each time 



The self-similarity structure of the rule 150 elemen- 
tary cellular automaton generated space-time fractal is 
qualitatively different from the Sierpinski triangle gen- 
erated by rule 90. While the iteration itself generalizes 
the concept of a Fibonacci iteration to vectors of grow- 
ing dimension, the blockwise sum exactly is given by the 
Fibonacci series multiplyed by a scaling factor 2". The 
iteration rule for the total activity derived here allows to 
compute the total activity without simulating the spatial 
dynamics, thus considerably eases the numerical compu- 
tation. 
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Hovever in contrary to the case of rule 90, generating the 
Sierpinski gasket, and where the cofficients of (1 + x) 
simply are the binomial coefficients, this does not lead 
to a closed expression for X(t) rule 150. 
The iteration in p| p. 614 gives wrong numbers for x( n )- 
Of order o(t 3 ), optimized o(t 2 ); the spatial simulation is 
also o(t 2 ); the iteration I14H is o(t) as the matrices are 
diagonal. 
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